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ABSTRACT 

Long-duration Gamma Ray Bursts (GRBs) are linked to the collapse of massive stars and their hosts 
are exclusively identified as active, star forming galaxies. Four long GRBs observed at high spectral 
resolution at redshift 1.5 < z < 4 allowed the determination of the elemental abundances for a set of 
different chemical elements. In this paper, for the first time, by means of detailed chemical evolution 
models taking into account also dust production, we attempt to constrain the star formation history 
of the host galaxies of these GRBs from the study of the measured chemical abundances measured in 
their ISM. We are also able to provide constraints on the age and on the dust content of GRB hosts. 
Our results support the hypothesis that long duration GRBs occur preferentially in low metallicity, star 
forming galaxies. We compare the specific star formation rate, namely the star formation rate per unit 
stellar mass, predicted for the hosts of these GRBs with observational values for GRB hosts distributed 
across a large redshift range. Our models predict a decrease of the specific star formation rate (SSFR) 
with redshift, consistent with the observed decrease of the comoving cosmic SFR density between z ^ 2 
and z — Q. On the other hand, observed GRB hosts seems to follow an opposite trend in the SSFR 
vs redshift plot, with an increase of the SSFR with decreasing redshift. Future SSFR determination in 
larger samples of GRB hosts will be important to understand whether this trend is real or due to some 
selection effect. Finally, we compare the SSFR of GRB050730 host with values derived with a sample of 
Quasar damped Lyman alpha systems. Our results indicate that the abundance pattern and the specific 
star formation rates of the host galaxies of these GRBs are basically compatible with the ones determined 
for a sample of Quasar damped Lyman alpha systems, suggesting similar chemical evolution paths. 

Subject headings: Gamma rays: burst. Galaxies: high-redshift. Galaxies: abundances; interstellar 
medium. 



1. INTRODUCTION 

Gamma-ray bursts (GRBs) are valuable tools to in- 
vestigate the high-redshift Universe. The extreme lumi- 
nosities of their afterglows have allowed detailed stud- 
ies of the physical properties of the interstellar medium 
(ISM) of their host galaxies. In particular, in a few cases, 
the determination of their chemical abundance pattern 
has b een possible (Savaglio ct al.' '2003"; 'Vree swiik et al. I 



2004 iProchaska et al. 2004; PenDrascct_aL 
Levesaue &: Kewlev I 120071 : IWiersema et al. 



2006; 



2007 



Prochaska et al. T 2007a ). as well as the determination 



of a few quantities related to their host galaxies, such 
as their dust con tent, their ste l lar m as s and their star 



form a tion rate (iBloom et aL _ 19981: Le Floc'h et al 



20021: iDiorgovski et al. I |2(p|; 



Le Floc'h et al. I 120061 : 
Savaglio et al.ll2008D . 



Chen et al. 



Christensen et al. 



ICastro Ceron et ~ 



2004; 



20061: 



"(|2005D first reported on the chemical abun- 
dances for the damped Lyman Alpha system (DLA) asso- 
ciated with the host galaxy of GRB 050730. Their analysis 
showed this gas was metal poor with a modest depletion 
level (see also Starhng et al. 2005). These results were 
subsequ ently expanded and tabulated bv lProchaska et al. I 
(|2007a^ (hereafter P07, see also D'Eha et al 2007) and we 
adopt th eir values adjusted to the solar relative abun- 
dances of iGrevesse et al. 1 ()2007[ ). 



In this paper, for the first time we aim at determining 
the star formation history of some GRB DLA systems 
by analysing the chemical abundances measured in their 
afterglow spectra. For this purpose, we use a detailed 
chemical evolution model, which allows us to predict the 
evolution of the abundances of the chemical elements stud- 
ied in the GRB DLAs of P07. Our aim is to constrain 
several relevant astrophysical quantities, such as the star 
formation rates, the dust content and the ages of the host 
galaxies of these GRBs. Our approach has already turned 
out to be fruitful to study Quasar DLAs on the basis of 
their chemical abundance pattern, allowing us to derive 
crucial info rmation on their nature and on their p hysical 
properties (jDessauges-Zavadskv et al. Il2004 I2007D . The 
plan of this paper is as follows. In Section 2, we briefly 
introduce the chemical evolution model. In Section 3, 
we present our results and in Section 4 we draw some 
conclusions. 

2. THE CHEMICAL EVOLUTION MODEL 

2.1. The dwarf irregular model 

The model used to derive the star formation histories 
of these GRB DLAs host galaxies i s similar to that devel- 
oped by iBradamante et alH (|1998f) to study dwarf irreg- 
ular galaxies. Since some chemical species studied in this 
paper are refractory (C, O, Si, Mg, Fe, Ni), in the chem- 
ical evolution model used here we include also a detailed 
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treatme nt of dust production and destruction following the 
work of ICalura et aFI (|2008D(CPM08, see section 2.2). 
We choose to use the model for dwarf irregular galax- 
ies since several observational investigations provided 
strong evidence that most of the GRBs originate in 
gas rich, star-forming sub-luminous {L < L^) galaxies 
with re latively low met a llicitie s {Z < Zir,) (iBloom et al. I 



19981: iFruchter et al. I 119991: iProchaska et al. I 
Vreeswiik et al. II2004D . 



20041: 



The dwarf galaxy is assumed to form by means of a contin- 
uous infall of pristine gas until a mass Mtot is accumulated. 
The evolution of dwarf irregular galaxies is characterized 
by a continuous star formation history. The star formation 
rate ip{t) is directly proportional the gas fraction G(t) at 
the time t, according to the Schmidt law ([Schmidt II1959I ) 

m = ^G(t). (1) 

where ly is the efficiency of star formation. We take into 
account the energy feedback from both Type la and Type 
II supernovae (SNe). They are responsible for the onset 
of a galactic wind, when the thermal energy of the ISM 
exceeds its binding energy, which is related to the pres- 
ence of a dark matter halo in which the galaxy is em- 
bedded (for more details, see Bradamante et al. 1998, 
Lanfranchi & Matteucci 2003). The thermal energy of 
the gas is controlled mainly by the thermalization efficien- 
cies of the supernovae explosions {rjsNii for SNe II and 
rjsNia for SNela) and stellar winds {risw)i which control 
the fraction of the energy restored to the ISM which is con- 
verted to thermal energy of the gas (see Bradamante et al. 
1998). The binding energy of the gas, on the other hand, 
is strongly influenced by assumptions concerning the pres- 
ence and distribution of dark matter (Matteucci 1992). A 
diffuse {Re/Rd—0-^, where Re is the effective radius of the 
galaxy and Rd is the radius of the dark matter core) but 
relatively massive (Mdark/MLum = 10) dark halo has been 
assumed for each galaxy. The time at which the wind de- 
velops depends on the assumed star formation efficiency. 
In general, the higher the SF efficiency, the earlier the wind 
develops. 

No instantaneous recycling approximation is adopted, i.e. 
the stellar lifetimes are properly taken into account. This 
allows us to compute in detail the contributions by low 
mass stars, Type la, Type II SNe to the chemical en- 
richment of the ISM. The main physical quantities related 
to these processes are the stellar yields, representing the 
fractions of processed matter restored by the stars into 
the ISM. The stellar yie lds used in this work are from 
iMevnet fc Maed eF* (^200 ^ for low and i nterm ediate mass 
stars (LIMS), from Wo oslev fc Weaverl (jl995D for massive 
stars and from Iw amoto et al. I (|1999D for Type la SNe. 
We assume a Salpeter initial mass function (IMF). With 
such an IMF, it is possible to reproduce the chemical abun- 
dances and the observation al features of dwarf irregulars 
(jCalura fc Matteucci 20061) and to accou nt for the local 
metal budget ( Calura fc Matteucci Il2004 ). We assume a 
cosmological model characterized by = 0.3, J^a = 0.7 
and a Hubble constant Hq = 70kms^^ Mpc~^ . 

2.2. N nucleosynthesis: indications from QSO DLAs 

Although N nucleosynthesis is still a matter of debate, 
most of N production arise from low and intermediate mass 
stars since massive stars produce only a small fraction of 



it (e.g. Chiappini, Matteucci fc Ballero 2005). In partic- 
ular, in the solar neighborhood the typical timescale for 
N enrichment should be 0.25-0.30 Gyr. However, the N 
produced in massive stars is important for understanding 
N abundances at very low metallicities. In principle, N 
produced in massive stars should be secondary, namely 
produced by means of the CNO cycle starting from the C 
and O abundances originally present in the star. In the 
last years Meynet fc Maeder (2000, 2002) and more re- 
cently Chiappini et al (2006) have shown that very metal 
poor stars can be fast rotators and that rotation can pro- 
duce primary N, namely N originating from the C and O 
produced by the star, not those originally present in the 
star at birth. In this CcLS6, cl high N abundance is expected 
at very low metallicities. It is not clear whether this high 
production of primary N in massive stars continue at in- 
termediate metallicities. In our case we adopt the stan- 
dard prescriptions for secondary N in massive stars and 
secondary plus primary N in low and intermediate mass 
stars. 

In Fig. [H we show the predicted evolution of log(N/0) 
vs time (upper panel) and vs metallicity (lower panel) 
assuming various star formation efficiencies. The mod- 
els are computed with the nucleosynthesis prescriptions 
used here (thick lines) and assuming primary N in mas- 
sive stars (thin lines), following Matteucci (1986). The 
predictions are compared with the observed N/0 values 
measured in high redshift QSO DLAs by various authors 
(Pettini et al. 2008 and references therein). The prescrip- 
tions adopted here are in good agreement with the most 
recent data published by Pettini et al. (2008) for DLAs , 
whose metallicities are not as low as metal poor halo stars 
where a high N abundance seems to be found (Spite et al. 
2005). The assumption of primary N in massive stars pre- 
dicts high N/0 values at low metallicity, i.e. immediately 
after the beginning of star formation. However, as can be 
seen from Fig.[TJ the highest N/0 values in DLAs can be 
explained also without assuming primary N production in 
massive stars, but by assuming a low star formation effi- 
ciency. 

Without primary N production in massive stars, the value 
log(N/0)=-2.6, corresponding to the lowest value ob- 
served in DLAs, is reached by our models after 0.05 Gyr 
with V — 0.01 Gyr~^ and after 0.1 Gyr with v = 1 Gyr"-'^. 
Finally, it is worth stressing that, as remarked also by Chi- 
appini et al. (2003), the stellar yields of Meynet fc Maeder 
(2002) do not include the third dredge-up and hot-bottom 
burning, thus their N yields for intermediate mass stars 
should be regarded as lower limits. 

2.3. Dust production and destruction 

The model for dust evolution used in this paper adopts 
the same prescriptions as CPM08. For a detailed descrip- 
tion of the dust evolution model, we address the reader to 
that paper. Here, we briefly summarize our assumptions. 
For the refractory chemical element labeled i, a fraction 
5f^, Sf"", and is incorporated in dust grains by low 
and intermediate mass stars, type la SNe, and type II SNe, 
respectively. These quantities are the dust condensation 
efficiencies of the element i in various stellar objects. Here 
we assume Sf^ = d("- = = Si = 0.1. This choice 
is supported by recent mid-infrared observations of one 



The star formation history of the host galaxies of GRBs 



3 



Supernova (Sugerman et al. 2006), which provided an up- 
per hmit of < 0.12. The value assumed here is also 
supported by theoretical studies of the local dust cycle 
(Edmunds 2001). Our assumption is compatible with val- 
ues provided by other theoretical and observational stud- 
ies (see Morgan & Edmunds 2003). However, other studies 
suggest that the dust condensation cfRciencies of SNe may 
be even lower, with values ranging from ^ 2 x 10^"' to 0.15 
(Zhukovska et al. 2008). 

CPM08 showed that the dust depletion pattern is strongly 
dependent on the star formation history. The depletion 
pattern observed in galaxies of different morphological 
types provides us with the most complete set of dust deple- 
tion data for a range of different star formation histories. 
The refractory elements involved in this study are C and 
O which in the local ISM show mild depletion levels, and 
Mg, Si , Ni, and Fe, which local ly show higher depletion 
levels (jSavage fc Sembacli1ll996[ ). A debate is currently 
open on the nature of S and on whether it may be in- 
corporated into dust grains or not. Observational studies 
of the depletion patterns of the cold ISM of the Milky 
Way report for S dep letion values between and -0.2 dex 
(jWeltv et al. Ill999al) . More recent data, based on observa- 
tional determinations of the gas-phase abundances of some 
chemical elements in the Local Interstellar Cloud, support 
a con siderable depletion for S of -0.3 dex (Ki mura et al. I 
|2003| ). However, these high depletion values have been de- 
rived by assuming as cosmic S abundance the solar value 
observed by Grevesse & Sauval (1998), i.e. As = 21.4 
per 10^ hydrogen atoms. The recent update by Grevesse 
et al. (2007) suggests a solar S abundance As = 14 per 
10^ Hydrogen Atoms, which implies a S dust fraction in 
the Local Interstellar Cloud lower than 0.4, i.e. a loga- 
rithmic depletion higher than -0.2 dex. However, other 
indirect arguments indicate that S may have a refractory 
nature. Detailed models of dark insterstellar clouds chem- 
istry suggest that sulphur must be depl eted by two to three 
orders of magnitude in dark regions (|Ruffle et al. I [19991 
IScappini et al. IF2003D . A recent investigation of dust and 
element abundances in planetary nebulae has shown that 
in these systems, de pletion is likely to be severe for various 
species, including S (jPhillips 1120071 ). In cometary dust par- 
ticles , S is present in iron sulphide (FeS) grains (jLodders I 
l2003f ). A broad feature ascribed to FeS emission has been 
recent ly seen in the infra red spectra of young stellar ob- 
jects ([Keller et al. ||2003 ). Furthermore, an indirect evi- 
dence of the fact that S may suffer depletion comes from 
the fact that its condensation temperature is considerably 
higher than the ones of C and O, which are known to be 
depleted in the gas phase (see Welty et al. 1999). From 
a photospheric abundance study of a RV Tauri star in the 
Large Magellanic Cloud, Reyniers & Van Winckel (2007) 
showed how the element abundance anti-correlates with 
the condensation temperature, possibly implying a corre- 
lation between the depletion and the condensation temper- 
ature. In the light of these facts, in our study we consider 
two different scenarios for dust grain composition. In the 
"expected dust" scenario, we assume that the only ele- 
ments going into dust grains are Ni and Fe. The reason 



of the name chosen for this scenario lies in the fact that 
mainly indirect evidences indicate that S could be a re- 
fractory element. We consider also an "alternative dust" 
scenario, where also S is refractory. 

Dust grains can be destroyed by the propagation of super- 
nova sho ck wa ves in the warm/io nised interstellar medium 
fMcKee i 19891: 1 Jones et anil994[ ). liGdust,t is the fraction 
of the element i locked into dust and G is the gas fraction, 
the destruction rate is calculated as Gdust.i/Tdestr, where 
Tdestr is the dust destruction timescale, calculated as: 

G 

Tdestr — (f^MsNR.)^ ^ ■ (2) 



(iMcKee I [19891: iDwek I I1998D . Msnr = 13OOM0 is the 
mass of the in t erstellar gas swept up by the SN remnant 
([McKee 1119891 iDwek et al. 1120071 ). Rsn is the total SNe 
rate, i.e. the sum of the rates of Type la and Type II SNe. 
Unless otherwise specified, we assume that no dust ac- 
cretion is taking place in the GRB host galaxies . This 
phenomenon occurs in j? 2-rich molecular clouds (|Dwek I 
ri998; Calura.et air[ l2008t) . Our choice is motivated by 
the fact that very little molecular H is observed in local 
dwarfs, with molecular-to-atqmic gas fractions of ~ 10% 
or lower ([Lisenfeld fc Ferrara |[T99alClavton et al. II1996D . 
Our assumption is further supported by the fact that a 
very little a mount of H2 observed i n the spectra o f GRB 
afterglows (IVreeswiik et al. I |2004 iFvnbo et al. I [20061: 



[Tumlinson et al. |[2007l: [Whalen et al. [[2008[ ). This might 
indicate that in the presence of intense SF, molecular 
clouds could rapidly dissolve, likely allowing very little 
dust accretion to occur. We will test the reliability of this 
assumption while discussing our results. 



3. RESULTS 

In this section, we attempt to constrain the star for- 
mation histories for the set of 4 GRB host galaxies. The 
abundance ratios measured by P07 for the 4 GRB hosts are 
presented in Table 1, as well as the predicted abundances 
computed by means of the best model, when possible. 
We address the reader to the paper by [Prochaska et al. 1 
f2007a) for further details on the GRI3 observations and 
on the techniques used to derive the ISM abundances of 
the GRB host galaxies. 

3.1. Abundance ratios as a diagnostic to infer the star 
formation history 

In chemical evolution models, the abundance ratios be- 
tween two elements formed on different timescales can be 
used as "cosmic clocks" and provide us with information 
on t he roles of LIMS and SNe in the enrichment pro- 
cess ([Matteucci Il2001[ ). In particular, the study of abun- 
dance ratios such as [a/Fe] and [N/a] [3 is quite useful, 
since the a-elements (O, S, Mg) are produced on short 
timescales by Type II SNe, whereas the Fe-peak elements 
and nitrogen are produced on long timescales by Type 
la SNe and low and intermediate-mass s tars, respe c tively . 
Calura et ahl (|2003D : [Dcssaugcs-Zavad skv et al. I (|2004 
20071) have shown that the simultaneous study of the abun- 
dance ratios between different elements as functions of a 

^ All the abundances between two different elements X and Y are expressed as [X/Y] = log{X/Y) — log(X/Y)Q, where (X/Y) and {X/Y)q 
ar e the ratios between the mass fractions of X and Y in the ISM and in the sun, respectively. We use the set of solar abundances as determined 
bv lGrevesse et al. I ll2007t) . 
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Fig. 1. — Upper panel : Predicted log(0/N) vs time for models adopting various star formation efficiencies. The thin lines include primary 
N in massive stars following Matteucci (1986). The thick lines do not include primary N production in massive stars. The shaded area 
represents the range of log(N/0) values observed by various authors (Pettini et al. 2008 and references therein) in a sample of QSO Damped 
Lyman Alpha systems. Lower panel: Predicted log(0/N) vs log(0/H)+12. Thick and thin lines as in the upper panel. The solid squares are 
observational values in high-redshift QSO DLAs (Pettini et al. 2008 and references therein). 



Table 1 



Ab. ratio 


GRB 050730 


GRB 050820 


GRB 051111 


GRB 060418 




Obs 


Fred 


Obs 


Fred 


Obs Fred 


Obs 


Fred 


[S/H] 


-2.26 


-2.34 


-0.63 


-0.76 






-1.27 


[Fe/H] 


-2.54 


-2.50 


-1.64 


-1.81 


-1.78 


-2.26 


-2.22 


[C/Fe] 


>-0.91 


-0.23 


>-0.34 


-0.42 






-0.32 


[N/S] 


-0.75 


-0.91 


>-0.57 


-0.79 






-0.88 


[S/Fe] 


+0.28 


0.16 


+1.01 


1.05 






0.93 


[Si/Fe] 


>-0.28 


0.15 


>+0.45 


-0.01 


>+0.72 


>+0.61 


0.09 


[0/Fe] 


>-0.91 


0.35 


>-0.32 


0.07 






0.34 


[Mg/Fe] 


<-l-0.90 


0.18 


+0.93 


-0.10 


>-0.98 


>-0.97 


0.03 


[Ni/Fe] 


-0.12 


-0.12 


+0.10 


-0.19 


-0.09 


-0.11 


-0.16 


[Zn/Fe] 




-0.05 


+0.95 


0.94 


>+1.14 


+0.65 


0.74 



Abundance ratios used in this work as measured in four GRB host galaxies and as predicted by means of our best models. All the 
abundances are normalized to the solar values, from Grevesse et al. (2007). 



metallicity tracer (such as [S/H] or [Fe/H] ) can be used a given system, whereas the study of the abundance ratios 
to constrain the nature and the star formation history of 
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Fig. 2. — Observed abundance ratios versus metallicity for the host galaxy of GRB 050730 as derived bv lProchaska et al. I l)2007af ) (empty 
squares with error bars and solid triangles). The thin solid line, dotted line and dashed line are predictions computed by means of a chemical 
evolution model for a dwarf irregular galaxy assuming for the star formation efficiency v = 0.005Gyr~^ , v = O.QlGyr~^ , and u = 0.02Gyr~^ , 
respectively, with "expected dust" composition. The thick lines are as above, but with "alternative dust" composition, i.e. with S depleted. 



1 I I -1 




2345672 3 4 5 
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Fig. 3. — Observed abundance ratios versus redshift for the host galaxy of GRB 050730 as derived bv IProchaska et al. I Il2007al) (empty 
squares with error bars and solid triangles). The thin solid line, dotted line and dashed line are predictions computed by means of the best 
model {u = 0.01Gyr~^ , "expected dust" composition) and assuming zj^ = 4, 2j = 5 and Zf = 6, respectively. 
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versus redshift can be used to derive constraints on the 
age of the system. 

In these 4 GRB DLAs, P07 have determined the abun- 
dance ratios for various elements, with typical errors of 
~0.1 dex. For several other elements, only upper or lower 
limits were derived for the abundance ratios (see Tab. 1). 
We use the complete set of elements studied by POT in 
order to constrain the SF history of the host galaxies of 
these 4 GRBs. In the following, we will present the results 
of our analysis, case by case. 

3.2. GRB 050730 

In this Section, our aim is to constrain the main pa- 
rameter of the star formation history of the host galaxy of 
GRB 050730, i.e. its star formation efhciency i'. We need 
to find the model which reproduces at best the observed 
abundance ratios. 

In this study, the only parameter allowed to vary is the 
star formation efficiency. Its value is set by requiring that, 
with the adopted SF efficiency, the observed abundance 
ratios are reproduced. 

We assume various values for the star formation effi- 
ciency ly, spanning from Vmin — O.OOSGyr^l to Vmax — 
O.lGj/r"!. In figure 2, we show the predicted evolution of 
various abundance ratios vs metallicity, traced either by 
[S/H] or by [Fe/H], calculated for a few models assuming 
different SF efficiency values and compared with the ob- 
served abundance ratios. 

Any model with a star formation efficiency v < 
0.005Gyr~^ is rejected, since they match the abundances 
of the non-refractory element N at ages larger than 1.53 
Gyr, which is the age of the Universe corresponding to 
the redshift of the GRB zgrb = 3.96. In Fig. 2, the 
thin lines represent the predictions computed with the "ex- 
pected dust" scenario, i.e. where S is undepleted whereas 
the thick lines are the predictions within the "alternative 
dust" scenario, where S is depleted. 

Some of the shown observational abundance ratios are 
lower limits ([C/Fe, 0/Fe] or upper limits [Mg/Fe], and do 
not provide any strong constraint which can help us decide 
which model is the best in reproducing the data. Useful 
information come from the [N/S] vs [S/H] and [S/Fe] vs 
[Fe/H] plots. The [N/S] vs [S/H] plot tells us that, with- 
out taking into account S depletion, no model allows us to 
reproduce the observed [N/S] value. On the other hand, a 
small S depletion of ~ 0.1 dex is enough to have the ob- 
served [N/S] reproduced by the two models with v = 0.005 
Gyr~^ and with v = 0.01 Gyr~^. However, as explained 
in Sect. 12.21 our predicted [N/S] should be regarded as a 
lower limit, since the yields of Meynet & Maeder (2002) for 
low and intermediate mass stars neglect the third dredge 
up, which may play an important role in N production 
(van den Hoek & Groenewegen 1997). 
The [S/Fe] vs [Fe/H] plot indicates that with S undepleted, 
the two models with i/ = 0.01 Gyr~^ and = 0.02 Gyr^^ 
reproduce the observed [S/Fe], however, once we take into 
account S depletion, only the model with i/ = 0.02 Gyr^^ 
is consistent with the data. On the other hand, the [Ni/Fe] 
vs [Fe/H] plot does not allow us to disentangle among the 
different models, since models with various SF efficiencies 

^ To compute the observed dust-to-gas and dust-to-metal ratios, wc assume for the Small Magellanic cloud a DgMC 
et al. 1990, Rubio et al. 2004) and Zq = 0.012 (Grevesse et al. 2007). 



present very similar [Ni/Fe] values. 

From the present study, we conclude that, within the "ex- 
pected dust" scenario, i.e. neglecting any S depletion, the 
best model is the one with ly — 0.01 Gyr~^, with the small 
discrepancy between our predictions and the models pos- 
sibly due to the fact that we are underestimating [N/H] 
since we assume no N production in massive stars. 
Within the "alternative dust" scenario, the best model is 
the one with i/ = 0.02 Gyr^^. Any model with SF effi- 
ciency ly = 0.03 Gyr~^ or higher leads us to underestimate 
the [N/S] value by ^ 0.2 dex. This can maybe reconciled 
by assuming primary N production in massive stars and 
more accurate yields for low and intermediate mass stars, 
including also the third dredge up. However, since a fur- 
ther investigation of the N is beyond the aim of the present 
paper, we think it is reasonable to reject any model under- 
estimating [N/S] by more than 0.2 dex , corresponding to 
the total error bar of the measured abundances considered 
here. 

Our conclusion is that, within the "expected dust sce- 
nario", the best model is the one with with ly = 0.01 
Gyr~^, whereas within the 'alternative dust" scenario, the 
best model is the one with i/ = 0.02. 

Once we have constrained the SF efficiency of the GRB 
050730 host galaxy, by studying its abundance ratios vs 
redshift we can derive constraints on its age. In the fol- 
lowing, we will use the two best models of the "expected 
dust scenario" and "alternative dust scenario" to derive a 
value for the age of the host galaxy of the GRB. 
In Fig. O we show the predicted redshift evolution of vari- 
ous abundance ratios, compared with the abundance ratios 
observed for GRB 050730. We show the predictions com- 
puted by means of the best model of the "expected dust " 
scenario by assuming three different values for the redshift 
of formation Zf, defined as the redshift at which the star 
formation turns on in the GRB host galaxy. The GRB 
host abundances are best reproduced by assuming z/ = 6. 
The difference between Zf and the redshift of the GRB 
host galaxy zqrb provides an estimate of the age of the 
system. The best model of the "expected dust" scenario 
points towards an age of ~ 0.6 Gyr. With the best model 
of the "alternative dust" scenario, a nearly identical value 
for the age of the GRB host can be derived. 
At this age, the predicted specific star formation rate is 4.7 
Gyr~l, both for the best model of the "expected dust" and 
"alternative dust" scenario. 

We have also a chance to test whether our assumption 
for dust evolution are appropriate, in particular the hy- 
pothesis of having neglected dust accretion. At the age of 
0.6 Gyr, with the best model of the "expected dust" pic- 
ture we predict a dust-to-gas ratio D = 1.2 x 10^^, which 
is in very good agreement with the observational upper 
hmilEI of Prochaska et al. (2007b) of 8 x 10"^. The pre- 
dicted dust-to-metals ratio is Dz — 0.03, again in agree- 
ment with the upper limit of 0.12 derived observationally 
by Prochaska et al. (2007b). Very similar values may be 
obtained by means of the best model of the "alternative 
dust" scenario. We conclude that the abundances of this 
system are not strongly influenced by dust depletion. Our 
assumptions concerning the treatment of dust evolution. 
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Fig. 4. — Observed abundance ratios versus metallicity for the host galaxy of GRB 050820 as derived bv lProchaska et al. ] ^M7^ (empty 
squares with error bars and solid triangles). The thin solid line, dotted line and dashed line are predictions computed by means of a chemical 
evolution model for a dwarf irregular galaxy assuming for the star formation efficiency u = 0.1Gyr~^ , v = lGyr~^, and u = 10Gj/r~^, 
respectively, with "expected dust" composition. The thick lines are as above, but with "alternative dust" composition, i.e. with S depleted. 




Fig. 5.— 



Symbols as in Fig 4. Thick and thin lines as in Fig 4, but including dust accretion. 
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Fig. 6. — Observed abundance ratios versus redshift for the host galaxy of GRB 050820 as derived bv iProchaska et al. I ll2007al ) (empty 
squares with error bars and solid triangles) . The dashed line, the dotted line and the solid line represent the predictions computed by means of 
the best model { i' = 0.lGyr~^ , dust accretion) within the "expected dust" scenario assuming Zf > 4, and Zf = 6, and Zf = 10, respectively. 



in particular having neglected dust accretion in molecular 
clouds, seems appropriate for the study of this GRB host. 

3.3. GRB 050820 

In fig. 4, we show various observed interstellar abun- 
dance ratios vs metallicity for the host galaxy of GRB 
050820. The observed abundance ratios are compared with 
predictions computed by means of three chemical evolu- 
tion models with SF efficiencies v = 0.1 Gyr~^, v — 1 
Gyr"-"^ and = 10 Gyr^^. In this case, we see how, with 
a model assuming only dust destruction and no dust ac- 
cretion, the observed abundance ratios do not allow us to 
determine which model is the most accurate in reproduc- 
ing the data. The predictions for the model with = 0.1 
Gyr"^ allows us to reproduce the [C/Fe], [0/Fe], [S/Zn] 
and are very close to the observed [N/S] lower limit. How- 
ever, no single model can reproduce the remaining set of 
abundance ratios shown in Fig. 4. In particular, the very 
high [S/Fe], [Mg/Fe] and [Zn/Fe] abundance ratios are not 
reproduced, even by assuming a star formation efficiency 
1^ — 10 Gyr~^, a value generally used to describe elliptical 
galaxies (Calura, Matteucci & Vladilo 2003). The same 
conclusion is valid taking into account S depletion (thick 
lines in Fig. 4). This GRB DLA presents peculiar abun- 
dance ratios, such as a high value for [Zn/Fe] which, at 
such a [Fe/H], in general is interpreted as the signature of 
dust depletion (see e.g. Vladilo 2004). To investigate the 
possibility that dust depletion may play a non negligible 
role in the ISM of the GRB 050820 host, we consider a 
chemical evolution model which takes into account dust 
accretion. To model dust accretion, the prescriptions are 
the same as described in CPM08. The dust accretion rate 
is Gdust,i/Taccr, whcrc Taccr is the dust accretion timescale 

Taccr = To^i/{l — fi) (3) 



where 

For the timescale tq i, typical values span from ^ 5 x 10'' 
yr, of the order of t he lifetime o f a typical molecular cloud, 
up to ~ 2 X 10* yr l|Dwek Ill998[ ). In this paper, we assume 
that the timescale ro,i is constant for all elements, with a 
value of 5 X 10^ yr (CPM08, Inoue 2003). 
In Fig. 5, we show our results for the predicted abun- 
dance ratios of the three models with i/ = 0.1 Gyr"-'^, 
I' — 1 Gyr^^ , and v — 10 Gyr~^ taking into account 
also dust accretion. By including accretion, a larger num- 
ber of abundance ratios is now reproduced by the model 
with v = 0.1, computed assuming no S depletion (thin 
lines in Fig. 5). In particular, the [S/Fe] and [Zn/Fe] ra- 
tios are now reproduced by the model with v — 0.1 Gyr~^. 
The very high [Si/Fe], [Mg/Fe] and [Ni/Fe] are not repro- 
duced by the same model, and they are likely to be due to 
differential depletion effects. Differential depletion is the 
most likely explanation for peculiar [Si/Fe] and [Mg/Fe] 
ratios observed also in some DLAs (Dessauges-Zavadsky 
et al. 2002, D07). In this paper, to limit the parameter 
space, we have assumed that all the chemical elements are 
incorporated into dust grains in the same proportions. A 
detailed modeling of differential depletion is beyond the 
aims of the present work. Future investigation on this 
topic, also in connection to dust depletion in DLAs, is re- 
quired to shed light on the chemical species interested by 
differential depiction. 

From Fig. 5, we see that, in presence of dust accretion, the 
predicted abundance ratios between a non-refractory ele- 
ment and a refractory one, such as [Zn/Fe], are the highest 
for the model with the lowest SF efficiency. This counter- 
intuitive result is due to the fact that, given the low con- 
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densation efficiencies assumed here, with high star forma- 
tion efficiencies {v > 1 Gyr~^), dust grains hardly survive 
in the ISM even if accretion is present. This result is due to 
the fact that a high star formation efficiency implies a high 
dust destruction rate, which the process of accretion is not 
able to compensate. This is consistent with the results of 
CPM08, who found that dust destruction by SNe is more 
efficient than dust production. However, C PM08 assume d 
the same dust condensation efficiencies as iDwek I (|1998[ ) , 
considerably higher than the ones used in the present pa- 
per (e.g., for Si and Fe, CPM08 assumed 5si,Fe = 0.8). In 
this case, also a galaxy with a SF efficiency of ^ 10 Gyr"-'^, 
such as an elliptical galaxy of luminous mass IQ^^Mq, can 
present a very high dust fraction, up to 0.9-1, in the period 
when star formation is active. 

It is also interesting to note that, in the "alternative dust" 
scenario, including S depletion and accretion, no model 
provides a satisfactory fit to the observed data. 
Regarding GRB 050820, we conclude that the model with 
SF efficiency v — 0.1 Gyr^^, "expected dust" composition 
and dust accretion is the best in reproducing the observed 
abundance ratios vs metallicity. We now use this model 
to infer the age of the GRB 050820 host galaxy. 
In Fig. 6, we show various abundance ratios vs redshift as 
measured by P07 for GRB 050820, compared to our pre- 
dictions computed by means of the best model and assum- 
ing three different redshifts of formation: Zf — A, zj = 6, 
and Zf — 10. It is clear that the majority of the abun- 
dance ratios, i.e. all the abundance ratios except [N/S], 
[Si/Fc], [Mg/Fe] and [Ni/Fc], aU underestimated for the 
reasons described above, may be reproduced by assum- 
ing Zf > 6. In this way, we are able to derive a lower 
limit of 1.5 Gyr for the age of the host galaxy of GRB 
050820. This allows us to derive an upper limit of 2.5 
Gyr~^ for the SSFR. Now, we investigate whether our as- 
sumptions concerning dust provide a realistic description 
of the dust content of this GRB DLA. For the GRB 050820 
DLA, Prochaska et al. (2007b) measure a dust-to-gas ra- 
tio 0.0008. Our results imply that, at the ages larger or 
equal to 1.5 Gyr, the predicted dust-to-gas ratio for the 
best model is > 0.00057, in very good agreement with the 
observed value. The predicted dust-to-metals ratio is 0.73. 
By assuming for the observed metallicity [Z/H] = [S/II]=- 
0.67 ±0.1 and by neglecting te uncertainty in the observed 
dust-to-gas, the measured dust-to-metals ratio is between 
0.23 and 0.36. Given the uncertainties of the parameters 
involved in our computation, in particular the dust con- 
densation efficiencies and the accretion timescale, which 
may even be larger by a factor of 2 than the value we have 
adopted, i.e. tq = 0.05 Gyr, we think our models provides 
a reasonable description of the dust properties of GRB 
050820. 

3.4. GRB 051111 

Of the elements studied in this paper, for GRB 051111, 
only 4 measured abundance ratios are available (see Tab. 
1). However, of these 4, 3 are lower limits and only one 
is a real measure. In fig. 7, we show the observed inter- 
stellar abundance ratios vs metallicity for the host galaxy 
of GRB 051111, compared with predictions computed by 
means of three chemical evolution models with SF efh- 
ciencies v = Q.l Gyr~^, v = 1 Gyr~^ and v = IQ Gyr^^. 



The predictions are shown both taking into account dust 
accretion (thick lines) and not (thin lines). The poor num- 
ber of abundance ratios measured for this GRB DLA does 
not allow us to constrain its star formation history. How- 
ever, the very large [Zn/Fe] is indicative of strong dust 
depletion. Unfortunately, no measures are available for 
the dust-to-gas and dust-to-metals ratios, hence it is not 
possible to assess accurately the amount of dust depletion 
characterizing this system. 

3.5. GRB O6O4I8 

Also For GRB 060418, only 4 measured abundance ra- 
tios are available to be compared with our predictions. 
In fig. 8, we show the observed interstellar abundance 
ratios vs metallicity for the host galaxy of GRB 060418, 
compared with predictions computed by means of three 
chemical evolution models with SF efficiencies v = 0.1 
Gyr^^, V ^ \ Gyr^^ and v ^ 10 Gyr^^, taking into ac- 
count dust accretion (thick lines) and not (thin lines). In 
this case, we have two real measures and two lower limits. 
Also this system presents a supersolar [Zn/Fe] ratio, which 
is likely indicating the presence of dust depletion. The two 
actual measures allow us to derive some constraints on the 
star formation history and on the age of the host galaxy of 
GRB 060418. From Fig. 8, we see that the model which 
best reproduces the observed abundance ratios is the one 
characterized by v = 0.1 Gyr~^ and including dust ac- 
cretion. With this model, it is possible to reproduce two 
measured abundance ratios, the lower limit for [Mg/Fe] 
is compatible with our predictions, whereas the [Si/Fe] is 
not, probably owing to the effects of differential depletion, 
as explained above. We use this model to derive a value 
for the age and for the SSFR of the host galaxy of this 
GRB. 

In Fig. 9, we show the measured abundance ratios vs 
redshift, compared to our predictions computed with the 
best model, assuming three different values for the red- 
shift of formation zf. Zf = 1.6, Zf = 1.8 and Zf = 2.. The 
observed abundance ratios are reproduced by assuming 
Zf — 1.8. The redshift of the GRB is zgrb =1.5, which, 
with the cosmology adopted in this paper, implies an age 
of 0.63 Gyr and a SSFR of 2.8 Gyr'^. At this age, we 
predict a dust-to-gas ratio of 0.0005 and a dust-to-metals 
ratio of 0.71. Unfortunately, no observational value for 
the dust-to-gas ratio is available for this system. 

A summary of the predicted physical properties of the 
GRB host galaxies studied in this paper is presented in 
Tab. 2. 

3.6. Specific star formation rates and implications for 
cosmic chemical evolution studies 

The specific star formation rate (SSFR) is defined as the 
ratio between the SFR and the stellar mass. Its determina- 
tion does not require any a-priori assumption on unknown 
quantities involved in our study, such as the mass or the 
physical dimension of the system, which on the other hand 
have to be assumed in order to provide an estimate of the 
SFR or the SFR surface density. 

The SSFR provides an indication of the intensity of the 
star formation of the host galaxy, whereas the inverse of 
the SSFR represents the star formation timescale, hence 
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Fig. 7. — Observed abundance ratios versus metallicity for the host galaxy of GRB 051111 as derived bv lProchaska et al. ] ll2007a[) (empty 
squares with error bars and solid triangles). The thin solid line, dotted line and dashed line are predictions computed by means of a chemical 
evolution model for a dwarf irregular galaxy assuming for the star formation efficiency u = O.lGyr^^ , u = IGyr^^, and u = 10Gj/r~^, 
respectively, with "expected dust" composition. The thick lines are as above, but taking into account dust accretion 
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Fig. 8. — Observed abundance ratios versus metallicity for the host galaxy of GRB 060418 as derived bv lProchaska et al. ] lIloOTi) (empty 
squares with error bars and solid triangles). The thin solid line, dotted line and dashed line are predictions computed by means of a chemical 
evolution model for a dwarf irregular galaxy assuming for the star formation efficiency u = 0.1Gyr~^ , v = lGyr~^, and u = 10Gyr~^, 
respectively, with "expected dust" composition. The thick lines are as above, but taking into account dust accretion 
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Fig. 9. — Observed abundance ratios versus redshift for the host galaxy of GRB 060418 as derived bv iProchaska et al. I ll2007al) (empty 
squares with error bars and soUd triangles). The solid line, the dotted line and the dashed line represent the predictions computed by means 
of the best model { u = 0.1Gyr~^ , dust accretion) within the "expected dust" scenario assuming Zf > 1.6, and Zf = 1.8, and zj = 2, 
respectively. 



a high value for the SSFR i s associated to a young galaxy 
(jCastro Ceron et al. I [2006h . It is interesting to com- 
pare the SSFR values predicted by means of our models 
with the SSFR observed for a set of GRB host g alaxie s 
at redshifts < z < 2.7 b v lCastro Ceron et al. I (|2008l ). 
ICastro Ceron et alTl (|2008l ) used the Spitzer images of 30 
GRB hosts and estimated their stellar masses and their 
unobscured star formation rates by means of their K-band 
fluxes and rest-frame UV spectra, respectively. By means 
of these data, it has been possible to determine the SSFRs 
for the host galaxies of their sample. In the left panel of 
Fig.[ini we show the SSFR as a function of the stellar mass 
for two models with SF efficiency i' — 0.01 Gyr~^ (solid 
lines) and u = 0.1 Gyr~^ (dashed lines). These models 
have turned out to be the best in describing the chemi- 
cal abundance patterns of the GRB DLAs studied in this 
paper. The models shown in Fig. 10 assume various val- 
ues for the mass normalization Mtot'- Mtot = W^Mq (thin 
lines) and Mtot = lO^^M© (thick lines). This stellar mass 
ra nge brackets the values obse rved in the bulk of the data 
bv lCastro Ceron et al. (I2008D. The models are compared 
to the data bv lCastro Ceron et al. I ^00§). 
Until the onset of the galactic wind, each pair of models 
with the same SF efficiency plotted in Fig. (TUj is charac- 
terized by the same chemical evolution and by the same 
abundance ratios. This is due to the fact that in chemical 
evolution models, to a first approximation, the abundance 
ratios do not depend on the total mass of the system, which 
plays a key role only in the presence of outflows. 
The predictions shown in Fig. 1 10! have been computed un- 
til the present time, corresponding to 13.5 Gyr for the 
cosmology adopted in the present paper. This is done to 
be consistent wi th the fact that the re d shift r ange spanned 
by the data by ICastro Ceron et al. I (|2008f ) is large, en- 



compassing ^ 82% of the cosmic time. Several observa- 
tional values plotted in Fig. [10] represent upper limits for 
the stellar mass and lower limits for the SSFR. In fact, 
ICastro Ceron et al. I ()2008[ ) showed that dust extinction 



must be present in at least 25% of the GRB hosts, causing 
the SFR determinations derived from UV/optical observa- 
tions to be regarded as low er limits. The SSFR values by 
ICastro Ceron et al~l (|2008t ) should be regarded as lower 
limits to the actual values also owing to a "dilution" effect 
present in mid-infrared photometry studies. This effect is 
related to the fact that some hosts of the sample are not 
spatially resolved, hence, the SSFRs are computed by di- 
viding the UV SFRs by the total stellar masses, instead of 
dividing the SFRs by the stellar masses of the star forming 
regions, which are likely to encompass a small part of the 
host galaxy. For these reasons, if the corrections for dust 
extinction and dilution were performed and if the stellar 
masses were determi ned with_higlicr precision, at least half 
of the data of lCastro Ceron et al. I (|2008f ) should move to- 
wards our predictions. 

In Fig. [ini right panel, we show the redshift evolution of 
the SSFR for o ur models, and for the GR B host galaxies 
as observed bv ICastro Ceron et al. I (|2008 ). In this case, 
in all our models, we have assumed that star formation be- 
gins at the same redshift z = 5. Some of the predictions for 
the various models coincide across a large redshift interval. 
This is due to the fact that the SSFR is not dependent on 
the SF efficiency v since, to a first approximation, both 
the SFR and the stellar mass are proportional to the SF 
efficiency v. In some cases, the predictions diverge owing 
to the onset of a galactic wind. This event concerns the 
models with the lowest mass normalizations, represented 
by the two thin lines in the right panel of Fig. [TOl charac- 
terized by the shallowest gravitational potentials. As soon 
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as the galactic wind develops, a considerable amount of 
gas becomes unavailable for star formation. The time the 
wind develops depends on the mass of the system and on 
the star formation efficiency v. In general, the higher is 
the adopted v and the lower is the total mass, the earlier 
the wind develops. For this reason, after the onset of the 
wind, the SSFR does depend on the adopted star forma- 
tion efficiency v. 

According to our predictions, the SSFRs decrease with 
decreasing redshift. Our result in not surprising, since 
any galactic evolution model predicts an increase of the 
stellar mass and a decrease of the SFR with decreas- 
ing redshift. This is in agreement also with the ob- 
served cosmological evolution of the SFR, in particular 
with the observed decline of the cosmic star formation 
history with redshift at z < 2 (Hopkinsfe Beacom 2006), 
likely due to progressive gas consumption i n late-type spi- 
ral discs (ICalura fc Matteucci I l2003l l2006l). Furthermore, 
iPapovich et al. ( 20061 ) showed that the integrated specific 
SFR, defined as the ratio between the comoving SFR den- 
sity and the comoving stellar mass density, is a decreasing 
function of reshift. These results are apparently in contrast 
with the observed redshift evolution of the SSFR in GRB 
host g alaxies. In fact, the data by ICastro Ceron et alTl 
(|2008[ ) suggest an increase of the SSFRs with decreasing 
redshift (see the solid straight line plotted in [T0| . at vari- 
ance with any results concerning the redshift evolution of 
star for ming galaxies at red s hifts < z < 3. As stated 
also bv lCastro Ceron et al. I (|2008t ). the GRB host galax- 
ies with the most vigorous star formation are likely to be 
heavily dust obscured, preventing the localization of their 
afterglows. This is one possible reason why, in the right 
panel of Fig. [101 no objects lay in the right upper corner 
of the plot. In the future, it will be important to extend 
the samples of GRB hosts with known star formation and 
stellar mass measures, in order to understand the reason 
for this apparent discrepancy and understand whether this 
trend is peculiar to GRB hosts or due to some selection 
effect. Surveys are ongoing to comprise GRB host galaxy 
samples based on localizations using the X-ray afterglow 
which is less susceptible to dust (Fynbo et al. 2007). 

3.7. A Comparison with SSFRs of Quasar DLAs 

In this section, we compare the SSFR of our mod- 
els for GRB hosts with the values determined by means 
of the ch emical evolution models fo r QSO DLAs pre- 
sented by iDessauges-Zavadskv et al."l ()2007[ ). We focus 
on the set of dwarf irregular galaxy models characterized 
by a continuous star formation history, which best re- 
produce the abundances o f the s et of DLAs studied by 
iDessauges-Zavadskv et al.1 (|2007[ ). By means of these 
models, the SSFRs are calculated at the time corre- 
sponding to the age predicted for the DLAs. In this 
way, for the 5 DLAs towards QSOs Q B0841-hl29, PKS 
1157-F014, Q B1210-hl75, Q B2230-h02, Q B2348-1444, 
we find SSFR between 0.24 Gyr"! and 7 Gyi-\ In Table 
3, we show the SSFR values com puted for the 5 DLAs of 
IDessauges-Zavadskv et alTl (|2007t ) , compared to the values 
obtained for the hosts of GR B 050730. The SSFR com- 
puted for the DLAs studied bv lDessauges-Zavadskv et al. I 
((200I are in substantial agreement with independent 
estimates for a larger sample of DLAs analysed by 



iHenrv fc Prochaska I (|2007[ ). 

Although the technique used bv lDessauges-Zavadskv et al. I 
(|2007f ) to correct the observations for dust depletion is dif- 
ferent than the ones used here, it is interesting to note 
that the SF efficiencies of the models used to describe 
the abundances of QSO DLAs range from 0.05Gyr~^ 
to QAGyr^^ (see Table 7 of Dessauges-Zavadsky et al. 
2007), i.e. basically similar to the values used here. A 
quantitative analysis of the differences in the dust de- 
pletion corrections comput e d her e and the ones used by 
IDessauges-Zavadskv et al.1 (|2007[ ) is beyond the scope of 
the present paper, but it is certainly an interesting subject 
for future studies. 

It is may be interesting to compare the SSFR values 
for the QSO DLAs of D07 with the values found in this 
paper for GRB DLAs. The abundances of the DLAs stud- 
ied by D07 required various depletion corrections. To 
correct for dust depletions, D07 used the technique de- 
scribed by Vladilo (2004). The implied depletion correc- 
tions (considering the [Fe/H] abundance) were 0.14 dex for 
QB0841+129 at z=2.375, 0.28 dex for PKS1157-)-014, 0.14 
dex for QB1210-I-175, 0.33 dex for QB2230-I-02 and 0.12 
dex for QB2348-1444. 

The approach to deal with dust depletion used here is dif- 
ferent than the one by D07. However, the two approaches 
are consistent, since both predict high [Zn/Fe] ratios for 
extremely dusty systems, and solar or nearly solar [Zn/Fe] 
for dust-free systems. Moreover, possible depletion effects 
involving S do not influence our results (see Sect. 13. 2p . 
It is probably wise and cautious to compare the SSFRs ob- 
tained for the three most dust poor DLAs of D07, with the 
star formation history derived here for GRB 050730, for 
which both the observations and our predictions indicate 
a very low dust content. To perform the comparison, we 
choose the SSFRs derived for QB0841-hl29 at z=2.375, 
QB1210-hl75 and QB 2348-1444. The values found for 
these 3 QSO DLAs of IDessauges-Zavadskv et al. I (j2007t) 
are plotted in Fig. [TT] As can be seen from Fig. [Til 
the SSFR of these 3 QSO DLAs are higher than the 
SSFRs of the best models of the "expected dust" sce- 
nario (within the "alternative dust" scenario, a very sim- 
ilar curve is obtained). The majority of QSO-DLA sight- 
lines do not penetrate highly SF galaxies, being the most 
metal rich, most vigorously star-forming systems very rare 
in QSO DLA samples, possibl y owing to dust obscura- 
tion of the backgrou nd quasars ( Hopkins fc Beacom II2006I : 
I Vladilo et al. I [20081) . Current expectation is that QSO- 
DLAs, which are drawn according to HI cross-section, 
most frequently correspond to sub-L, galaxies, as evi- 
denced by their significantly sub-solar metallicities (e.g. 
Fynbo et al. 2008). At high z, the luminosity function 
is sufficiently s teep that sub-L.^ ga laxies dominate the in- 
tegrated SFR (jFvnbo et al. II2008D . For this reason, DLA 
samples should be able to trace the bulk of star formation. 
iProchaska et al. I (|2007af) discuss how GRB DLAs prefer- 
entially probe gas associated to the innermost regions of 
galaxies. This gas is generally denser and more metal rich 
than the gas probed by QSO DLAs which, on the other 
hand, are likley to arise in the outer galactic regions. If 
GRB DLAs probe regions of the Universe different than 
the ones of QSO DLAs, the GRB DLA and QSO DLA 
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Table 2 

Summary of the physical properties of the 4 GRB host galaxies studied in the present paper 



GRB 
050730 
050820 
051111 
060418 



Age (Gyr) 
0.61 
>1.53 

0.66 



SSFR {Gyr-') 
4.7 

< 2.5 Gyr-i 



Dust-to-gas 
1.2 X 10-6 
> 0.00057 

0.0005 



Dust-to-metals 
0.03 
0.73 

0.71 



Table 3 

Specific Star Formation Rates predicted for GRB hosts and QSO DLAs 
DLA SSFR {Gyr-') 



QB0841-I-129 


1.5 


PKS1157+014 


0.57 


QB1210+175 


0.48 


QB2230-I-02 


0.24 


QB2348-1444 
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GRB 


SSFR {Gyr-') 


GRB 050730 


4.7 


GRB 050820 


< 2.5 


GRB 050820 


< 2.8 
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Fig. 10. — (a): SSFR as a function of the stellar mass for two models with SF efficiency v = 0.01 Gyr~^ (solid lines) and v = 0.1 Gyr~^ 
(dashed lines), and assuming Mtot = 10^ Mp, (thin l ines) and Mtot = IO'^^Mq (thick lines). The solid squares and the triangles are the 
observational determinations bv lCastro Ceron et al. ] ll2008t) . The triangles mark lower or upper limits f or the derived stellar mass, SSFR or 
both, (b): Redshift evolution of the SSFR for our models, and for the GRB host galaxies as observed bv lCastro Ceron et al. IllHoOS). All the 
curves, squares and triangles are as in panel a. The thick dot-dashed line represents the best fit to the observed values, computed by means 
of a straight line y = ax + b. All the other curves, squares and triangles are as in panel a. 



samples may provide complementary information on cos- 
mic chemical evolution. 

4. CONCLUSIONS 

In this paper, for the first time we have used the ele- 
mental abundances observed in the host galaxy of 4 long 



GRBs to constrain their star formation history and its age. 
Our method is based on the simultaneous study of the 
abundance ratios between various elements synthesized by 
stars on different timescales, by comparing the measured 
abundance ratios to the predictions of a detailed chemi- 
cal evolution model. From the simultaneous study of the 
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0.5 



Fig. 11. — Predicted evolution of the SSFR vs [Fe /H] for the best model for GRB 50730 host galaxy. The open squares are the values 
determined for the three QSO DLAs of the sample bv lDessauges-Zavadskv et al. I ||2007| ) presenting the lowest dust depletions. 



abundance ratios between various elements vs metallicity 
(traced by [Fe/H], [S/H] or [Zn/H]), we have been able to 
constrain the star formation history of some of the host 
galaxies. By studying the abundance ratios vs redshift, 
we have been able to constrain the redshift of formation 
and the age of the GRB host galaxies. We tested mod- 
els both with and without dust accretion. Some peculiar 
abundance ratios, such as very high [Zn/Fe] ratios, seem 
to indicate severe dust depletion and, to be explained, re- 
quire the inclusion of dust accretion in the models. Our 
main results can be summarized as follows: 

1) Our results point towards star formation efficiencies be- 
tween 0.01Gyr~^ and O.lGyr^^, ages for the system be- 
tween 0.6 and >1.5 Gyr. The observed dust content of 
some GRB DLAs is also satisfactorily accounted for. 

2) The predicted specific star formation rate values are 
between 2.5 Gyr~^ and 4.7 Gyr~^ . These values are com- 
patible with o bservational estimates from a very recent 
compilation bv ICastro Ceron et al. I (|2008f ). Our models 
predict a decrease of the SSFR with redshift, consistent 
with the observed decrease of the comoving cosmic SFR 
density between z ~ 2 and z = 0, in agreement with nu- 
merous observations of star forming galaxies at low and 
high redshift. On the other hand, apparently the detected 
GRB hosts follow an opposite trend in the SSFR vs red- 
shift plot, with a slight increase of the SSFR with de- 
creasing redshift. In the future, it will be important to 
asses whether this apparent trend is real or if it is due to 
some selection effect affecting the observation of GRB host 
galaxies. 

3) The study of the abundance pattern of the host galaxy 
of GRB 050730 indicates that the predicted specific star 
formation rate is compatible with the values found for 
the the set of DLAs studied bv lDessauges-Zavadskv et al. I 



(|2007t) . This implies that, the host galaxy of GRB 050730 
may have followed a chemical evolution path similar to the 
ones typical of QSO DLAs. 

Our study favors the hypothesis that long duration GRBs 
occur preferentially in low metallicity, star forming galax- 
ies, characterized by low star formation efficiencies {u < 
0.1 Gyr-i). 

The fact that GRBs occur preferentially in low metallicity 
systems is indicated also by the luminosity-metallicity rela- 
tion ( Wolf fc Podsiadl owski . 200 7) and by some progenitor 
models for GRBs, where low metallicity hampers large loss 
of angular momentum and mass and the progenitor can re- 
tain a rapid rotating core, as req uired in particular in the 
collapsar model for long GRBs (|MacFadven &: Wooslev I 
,19991 

In the future, in order to have a better understanding of 
the chemical evolution patterns observed in high-redshift 
galaxies, it will be important to refine the study of dust 
depletion in local and distant galaxies, in order to clarify 
whether elements such as S have a refractory nature or 
not. 

Furthermore, it will be important to extend the study car- 
ried on in this paper to a larger set of GRB DLAs, in order 
to have an insight into regions of the universe complemen- 
tary to those probed by QSO DLAs and by other types of 
astrophysical objects, detected by means of difi^erent tech- 
niques. 
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